Dynamical Correlations among Vicious Random Walkers 
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Nonintersecting motion of Brownian particles in one dimension is studied. The system is con- 
structed as the diffusion scaling limit of Fisher's vicious random walk. N particles start from the 
origin at time t = and then undergo mutually avoiding Brownian motion until a finite time t = T. 
In the short time limit t<T, the particle distribution is asymptotically described by Gaussian Uni- 
tary Ensemble (GUE) of random matrices. At the end time t = T, it is identical to that of Gaussian 
Orthogonal Ensemble (GOE) . The Brownian motion is generally described by the dynamical corre- 
lations among particles at many times t\,ti, • ■ ■ , £m between t = and t = T. We show that the 
most general dynamical correlations among arbitrary number of particles at arbitrary number of 
times are written in the forms of quaternion determinants. Asymptotic forms of the correlations in 
the limit N — ► oo are evaluated and a discontinuous transition of the universality class from GUE 
to GOE is observed. 

PACS number(s): 05.40.-a, 05.50.+q, 02.50.Ey 

The vicious walk model in which many walkers ran- 
domly move without intersecting with others was in- 
troduced by Fisher and applied to wetting and melt- 
ing phenomena Q. Recently it has attained a re- 
newed interest, since intimate relations to other research 
fields, such as the theory of Young tableaux in com- 
binatorics asymmetric exclusion process (ASEP) 
in noncquilibrium statistical mechanics [5|, polynuclear 
growth (PNG) model of surface physics pM| and the the- 
ory of random matrices j^,^] , have been revealed one after 
another and brought progress in the study of these topics. 

One of the reigning concepts of these new applications, 
the universality class, comes from the theory of random 
matrices. Gaussian Orthogonal Ensemble (GOE) and 
Gaussian Unitary Ensemble (GUE) universality classes 
originate in the real symmetric and complex hermitian 
structures of random matrices. They also appear in the 
new applications in quite natural ways, despite the ab- 
sence of underlying matrix structure. In the theory of 
Young tableaux, number permutation and involution cor- 
respond to GUE and GOE, respectively @. In ASEP, 
the initial and boundary conditions demarcate the uni- 
versality classes of the current fluctuation pd|,^2[ . The 
height fluctuation of surface growth on a flat substrate 
belongs to the GOE class and a droplet growth is de- 
scribed by GUE §. 

How the universality classes appear in the vicious walk 
model? The GOE universality class is realized when 
walkers take t steps under the nonintersecting condition. 
One possibility to observe the GUE class is to impose 
an additional condition that the nonintersecting walkers 
come back to the original position after 2t steps ||[5],[ll] . 
In analyzing the diffusion scaling limit of the vicious walk 



model, Katori and Tanemura |Uj,|14| recently noticed an- 
other possibility. The walker distribution depends not 
only on the observation time (the number of steps) t but 
also on the time interval T in which the nonintersect- 
ing condition is imposed. Suppose that all the walkers 
start from the origin at time t = 0. While the ratio of 
the time t and the nonintersection time interval T, t/T, is 
small, the walker distribution is asymptotically described 
by GUE. When the end of the nonintersection time inter- 
val arrives, namely at t/T = 1, the distribution becomes 
identical to that of GOE. This means that, while walkers 
randomly move in the time interval between t = and 
t = T, a transition from GUE to GOE takes place. 

In this Letter we analyze the dynamical correlations 
among the vicious walkers in the transition region. Uti- 
lizing the equivalence to a multimatrix model in quan- 
tum field theory, we evaluate the dynamical correlation 
functions among arbitrary number of walkers at arbitrary 
number of times. The asymptotic limit of the large num- 
ber of walkers will be evaluated and we will find that the 
transition becomes discontinuous. 

Let us consider N independent symmetric simple ran- 
dom walks on Z = {• ■ • , —2, — 1, 0, 1, 2, • • •} started from 
N distinct positions 2s\ < 2s 2 < ■ ■ ■ < 2s^, Sj S Z. 
The position of the j-th random walker at time k > is 
denoted by R s k 3 and we impose the nonintersecting con- 
dition 



<K 2 



< 



1 < Vfc < K. 



(1) 



If a possible random walk satisfies the condition (1), then 
it is called a vicious walk. Let V(R S ^ = ey) be the realiza- 
tion probability of the vicious walks, in which the N walk- 
ers arrive at the positions 2e\ < 2e 2 < ■ ■ ■ < 2ejv, ej G Z, 
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at time K. 

A simplification is obtained in the diffusion scaling 
limit: we set K — Lt,Sj = s/Lxj/2, e,- = \/~Lijj/2 and 
take the limit L — > oo. In Ref. fl3|| , it is clarified that 

hm i ^ 0O (VI/2) Jv ^(i?^ /2 = yi % /2) = f(t;{y 3 } \ 
{xj}), where 



Here (t M +i = T) 



/(*; {%} I I x j}) 



■det[e-^- yk)2/2t ]^ k=1 ,..^ 



N- 



(2) 



This function shows how the nonintersecting probabil- 
ity of the Brownian particles on the rescaled lattice 
Z/(\/Z/2) up to time t depends on the initial and final 
positions. Therefore, for the nonintersecting Brownian 
motions in the time interval [0, T], the transition density 
from the configuration {xj} at time s to {yj} at time t 
is given by Q 



<p T {s,{xj};t,{ yj }) = 



f(t ~ s; {yj} | {xj})Af(T - t; {yj}) 



M(T-s;{xj}) 



(3) 



where N(t; { Xj }) = f yi< ...< VN U j= i {%} | {xj}). 

Note that there is an obvious temporal inhomogcneity 
since RHS depends not only t — s but also T — s and 
T — t. 

The dynamics of the Brownian particles is described 
by the dynamical correlation functions. Let us denote 
the positions of the Brownian particles at a time tj by 
x\, x°2, ■ ■ ■ , x J N . Then the dynamical correlation functions 
among the particles at times t%, t-2, • • • , £m are defined as 



nlr 1 ■■■ t 1 • • • • • r M ■■■ r M ) 

r N N N 



dxf 



j=X " j=n 1 + l j=n M +l 

M-l 

x Po ({x°j}) n {^l;Wi, K n+1 »- (4) 

m— 

Here Po^x®}) is the initial distribution at to — 0. Let 
us now suppose that all the particles start at the origin. 
Namely, we set po({x®}) = Y\j S(x { j) and obtain 

p(x 1 , • • • , x. ni ; • • • ; x 1 , • • • , x nM ) 



3C 



J ii dx)...j n </n< +i 

W AT 



M 

X TTdct[ 3 m (x™,X™ +1 )],, fc=1 ,.., A r. (5) 



^x 2 /(2t 1 ) -(x-y) 2 /(2(t 2 -t 1 )) 



e -(x-yf/(2(t m+1 -t m )) 
g m (x, y) = 2 < m < M. 



m+1 ^ro 



(G) 



m— 1 



At this stage we notice a direct correspondence be- 
tween our problem and the multimatrix models in quan- 
tum field theory. Itzykson and Zuber firstly analyzed a 
two matrix model in which two hermitian random matri- 
ces were combined fll|[l6| . Mehta and Pandey then cou- 
pled a real-symmetric and a hermitian random matrices 
to devise a (1 + 1) matrix model, as a mathematical inter- 
polation of GOE and GUE After [t had been real " 
ized that multimatrix models were useful in the quantum 
field theory on random surfaces [Q, Eynard and Mehta 
invented a method to evaluate the correlation functions 
among the eigenvalues of combined M hermitian matri- 
ces po[ . As a further generalization, Nagao proposed an 
(M+1) matrix model in which one real symmetric and M 
hermitian matrices were combined and showed that the 
correlation functions were generally written in the forms 
of quaternion determinants [^l| . We can readily see that 
the above dynamical correlation functions among vicious 
walkers have the same forms as the eigenvalue correlation 
functions of the (M + 1) matrix model. 

For simplicity we set N even and summarize the 
quaternion determinant formulas in the following. In 
terms of (x m = x, x„ = y) 

G mn (x, y) 

5(x — y), m = n, 
g m (x,y), m = n-l, 

n-l n— 1 (7) 

H dXj JJ g l (xi,xi + i), m<n-l, 

j=m+l l—m 

we firstly define 

F mn (x,y) 

oo pz' 

dz 1 / dz{G m M+1 (x, z)G n M+1 (y, z') 

■oo J —oo 

- G n M+1 {y 7 z)G m M+1 (x,z')}. (8) 
Let us introduce an antisymmetric inner product 

= ljdxj dyF n (y, x)[f(y)g(x) - f(x)g(y)\ (9) 

and construct monic polynomials R\{x) = x k + • • • of 
degrees k so that they satisfy the skew orthogonality re- 
lations: 

(Rlj(x),B^ l+1 (y)) = -{Rl l+1 (x),Rlj(y)) = rjSji, 
{Rl^xlRUy)) = 0, (R\j +1 (x),R\ l+1 {y)) = 0- (10) 
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We then define functions R^ix) and &™(x), m 
2,3, • • -,M+ 1, as 



Rf(x) = J dyRl(y)G lm (y,x), 

$I"(x)= [ dyRT(y)F mm (y,x). (11) 



Now matrices D mn , I mn and S mn are introduced as 

(AT/2)-l 



D T = E 



^2fe (^J 1 ) ^2fe+l ( X P ) -^2fe+l (^j" ) -^2fc (^P ) 



/c=0 



jmn 



and 



(iV/2)-l 

E 



(AT/2)-l 



$- fc+1 (^)^ fc (xr) - <^(zr)<^ fc+1 (xr 



fc=0 

Further we define 



''A' 



S"} n , m>n, 



and 



j7 I SIT ~ Gmn ( x ?, x i), m < 



jmn jmn , jrimn ( rn n\ 

1 jl ~ 1 il + r \ X j J X l )■ 



For a self dual N x N quaternion matrix Q — [qjk], 
a determinant TdetQ, originally introduced into random 
matrix theory by Dyson is defined as 



Tdet Q = ^2(-l) N ~ l 11 tr(q ab q bc ■ 



■qda) 



(12) 



where P denotes any permutation of the indices 
(1,2,- • • , N) consisting of I exclusive cycles of the form 
(a^b^c~^---^d^a) and (— l)^ - ' is the parity 
of P. Note that tiq is equal to a half of the trace of the 
2x2 matrix representation of q. 

Let us suppose that the elements of quaternion ma- 
trices B^, fi,u = 1, 2, • • • , M, have the following 2x2 
representations : 



jl Jl 



3,1 = 1,2,- 



(13) 



Then Nagao's result |EM asserts that 

p( x i, ' ' ' ? x ni > ' ' ' > X l ' ' ' ' ' X n M ^ 
cx Tdet [B^in^nv)], n,u = l,2,< 



(14) 



where each block B^^n^, n v ) is obtained by removing 
the n^+l, • • • , N-th rows and n u +l, n„ + 2, • ■ • , N- 

th columns from B^ v '. 



Let us examine the consequences of the quaternion de- 
terminant formula. Now we remark that the skew or- 
thogonal polynomials R\ (x) are explicitly written as 



Ri{x)=e /2 j^c ikj H j (£)c 



■3/2 



(15) 



3=0 

where £ = t 1 /(2T-t 1 ), 



Cn = y/t n (2T-t n )/T, 



OL2k j 
&2k+l j = 2 



-2fc„2fcr 



-2A:- 



-1 2fc+l 



2fc+l j 



4fc(5 2 fc- 



1 3, 



and Hj(x) are the Hermite polynomials. 

We begin with the simplest case M = 1 and m = 1. 
Putting the explicit formula of i?^ (x) into the quaternion 
determinant expression and utilizing the asymptotic for- 
mula for the Hermite polynomials, we can readily obtain 



p(xl)oc J2N-(x\/ Cl ) 2 , \x{\<s/2N Cl (16) 

7TCl v 

in the limit N — > oo. Thus the walker density always has 
a semicircle shape (Wigner's semicircle law), while the 
width of the semicircle is dependent on time and scaled 
by c\. This result suggests that, after introducing a new 
rescaled variables A™ = x™ /c m , the vicious walk model 
in the diffusion scaling limit becomes equivalent to the 
matrix Brownian motion model (23-27 which is normal- 
ized so that the width of the semicircle is a constant. 

The asymptotic forms of the dynamical correlation 
functions of matrix Brownian motion models were eval- 
uated by Forrester, Nagao and Honner |p8| . In par- 
ticular, in the edge region of the semicircle, dynamical 
asymptotic correlations are described by the Airy func- 
tion Ai(x). We can directly reinterpret Forrester, Nagao 
and Honner's result in the context of vicious walk model. 
Let us introduce rescaled temporal and spatial variables 
T m and XT as 



t m = 1 - 



2N- 



T, 



2 l/2jVl/6 



(17) 



and take the limit N — > oo. With an appropriate nor- 
malization, the dynamical correlation functions asymp- 
totically have the same forms as in eq. (|l4| ) except that 
the quaternion matrices B^ are replaced by 



CM" 

b 3l 



3,1 = 1,2,- 



,N, 



(18) 



where 
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-i r poo j 

V% = 4 \J d^'Ai^ + S )-{e-^Ai(Xf + S )} 
-J dse- T » a Ai(X? + s)— {e- T - s Ai(X^ + s)} , 

pOO pOO 

Iff = J dse- T " s Ai{X^ -a) J dve-^ v Ai(X^ - v) 

pOO poo 

- / dse-^ Ai(X!f - s) / dve~ T " v Ai(X? - v) 

JO J s 



and 



with 



c/ty _ / S ji> M > v > 



poo 

57= / dse^- T »> M{Xf + s)Ai(Xf + s) 
Jo 

+ iAi(Xf) d S e-^ s Ai(X; - fl ), 



dse (7 



->Ai(Xf + s)Ai(Xf + s). 



From the above asymptotic result we can extract infor- 
mation on limiting behavior. Let us firstly take the limit 
t ^ — ► 00 with the time differences r M — r v fixed. It can 
be readily seen that this limiting procedure is equivalent 



to set Vf" 



T;, and the second term of S^, zeros. Then 



the quaternion determinant is reduced to an ordinary de- 
terminant and the asymptotic correlation functions be- 
come temporally homogeneous. We find that they are 
the dynamical correlation functions within the universal- 
ity class of GUE. The equal time correlation functions are 
described by the Airy kernel j2j],|3(J in the temporally ho- 
mogeneous region. Because of the rescaling (|l7|), it can 
be seen that the GUE universality class survives until 
time t very close to T: only when T — t ~ 0(N~ 1 ^ 3 ), 
the transition to GOE class occurs. Therefore we can 
conclude that the transition from GUE to GOE class is 
discontinuous in the limit N — > 00. 
The second interesting case is Xj 1 



-00 with the po- 



sition differences Xj — X'( fixed. This reproduces the 
asymptotic dynamical correlation functions in the bulk 
region. They are spatially homogeneous and the equal 
time correlations are equivalent to Pandey and Mehta's 
asymptotic result [M. In this bulk region we can also 
observe the discontinuous transition from GUE to GOE. 

In summary, nonintersecting Brownian motion of N 
particles in finite time interval < t < T was studied in 
one dimension. Dynamical correlation functions among 
many particles at many times were written in the forms 
of quaternion determinants. The asymptotic forms of the 
dynamical correlation functions in the limit N — > 00 were 
evaluated and the universality class transition from GUE 
to GOE was found to be discontinuous. The compact 



asymptotic formula was derived above only after taking 
the diffusion scaling limit in which much information con- 
tained in the original lattice model was suppressed. In 
order to fully understand the lattice vicious walk model, 
we need to study other scaling limits as well. Further 
studies in this direction should be made in future works. 
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